Experimental determination of mb{Mz) in DELPHI 
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The running mass of the b quark as defined in the MS renormalization scheme, m^, was measured at the Mz 
scale using 2.8 million hadronic Z'^ decays collected by the DELPHI experiment at LEP. The result is 

mt{Mz) = 2.67 ± 0.25(stat.) ± 0.34(frag.) ± 0.27(theo.)GeV/c^ 

which differs from that obtained at the T scale, by mi,{Mr /2)~mi,{AIz) ~ 1.49±0.52 GeV/c^ . This measurement, 
performed far from the bb production threshold, provides the first experimental observation of the running of the 
quark masses. 
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1. Introduction 

The quark masses are fundamental parameters 
of the QCD lagrangian which are not predicted 
by the theory. Their definition is not unique 
since quarks are not observed as free particles, 
but confined inside hadrons. Therefore, the use 
of dynamical expressions is mandatory in order 
to determine the masses of the quarks. The per- 
turbative pole mass, Mq, and the running mass 
in the MS scheme, mg, are two definitions among 
the most frequently used. 

In the renormalization procedure of QCD, once 
the divergent part is subtracted, we deal with fi- 
nite quantities for the strong coupling constant. 
Us, and the quark masses, rriq. However, these 
quantities are not fixed parameters, but their ac- 
tual value depends on the energy scale of the pro- 
cess being considered. In other words, and mi, 
exhibit a running property. 

The running of as has been verified many times 
by different experiments (for a review see [ |], || ) . 
Although the running of the quark masses is as 
basic to QCD as that of as, it has never been 
experimentally tested before. The empirical de- 
termination of the quark masses at different scales 
must be regarded as a fundamental test of QCD. 

Nowadays, the LEP I data allow measuring 
the b quark mass at the Mz scale (far away 
from the bb production threshold), thus testing 
its running. According to the nib evolution, dic- 



tated by the Renormalization Group Equation 
(RGE), a change on to;, of ~ 1.3 GeV is expected 
when varying the scale from Mx /2 up to Mz , as 
TO6(Mt/2) « 4.16 ±0.18 GeV/c2 [|. 

In e+e^ ^hadrons, the three-jet rate is an ob- 
servable sensitive to the b mass. When testing 
the universality of as by using the three-jet rate, 
it was observed that due to the mass effects, b 
quarks radiate a 3 to 5% less gluons than light 
quarks do [ ^. So, reversing the argument: as- 
suming the as universality allows measuring the 
b quark mass by studying the three-jet rate. 

The former Leading Order (LO) calculations 
of the three-jet cross section in e^e^ including 
mass terms [ |^, ^ are not appropriate in order to 
evaluate the b quark mass, as these calculations 
can not discern between Mb and nib- Recently, 
Next to Leading Order (NLO) expressions for the 
three-jet rate in e+e~ are available [0, ^, ^, 
thus, enabling the measurement of to^. 

The observable used in this work is: 



'(ye)/r,V 



where F 
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3, "'iv^)/r?°r" 
nitJ.) {bo{yc,n) + ^^biiy^n) 



and Ff„t~"«« denote the differen- 



tial three-jet and the total cross sections for b and 
light quarks (£ = u,d,s), rb{n) = {'mb{^) / M zY ■ 
The coefficients &o and bi are calculated at NLO 
in reference [ 0|. j/c is the jet resolution parameter 
and ^ the characteristic energy scale. 



2. Data analysis 

The analysis reported was performed over a 
sample consisting of 2.8 million Z° hadronic de- 
cays recorded by the DELPHI detector at LEP 
through the years 1992 to 1994. 

2.1. Event selection 

Those hadroni c ev ents that satisfied the flavour 
tagging (section |2.3|) and the jet reconstruction 
criteria were retained for the analysis. 

The jets of the event were reconstructed by 
means of the Durham jet finding algorithm [ [l0|. 
By choosing the appropriate ?/c, it was possible 
to force the reconstruction of just three jets in 
every event. Then, a set of quality cuts were 
applied over each jet (minimum charged multi- 
plicity, enough visible energy, jet direction in the 
barrel part of DELPHI and overall planar config- 
uration). These conditions had to be simultane- 
ously fulfilled by all three jets. 

A total of 1,074,860 and 294,509 events entered 
in the t and h sample respectively. 

2.2. Flavour definition 

The true flavour of the event was defined as 
that of the quark coupled to the Z in the Z — > 
gg vertex, without prejudice to the new flavour 
production that possibly occurred in the splitting 
of a bremsstrahlung gluon (e.g. g — > 66). This is 
the same criterion adopted in [ |^ , thus it enables 
a direct comparison of our result with that NLO 
calculation. 

2.3. Flavour tagging technique 

The signed impact parameter of all charged 
particles in the event was used to set up an al- 
gorithm which enables the flavour tagging [ |ll) . 
The significance of each track was obtained by 
weighing its impact parameter with its associ- 
ated error. Thus, a function, V ., was constructed 
in order to estimate the probability of having all 
particles compatible with being generated in the 
events' Interaction Point (IP). By construction, 
light quark events had an uniform distribution of 
V . In the 55 events, the decay of long lived B 
hadrons led to particles generated in secondary 
vertices far away from the IP, besides, biasing V 
towards low values. Accordingly, 66 events were 



selected by requiring V < 5-10^^ and ££ with 
V > 0.2. The purity attained for the 5 sample is 
approximately an 85% and its efhciency a 55%, 
while for £ quarks both are about the 80%. The 
contents on cc events in each tagged sample (6 
and £) were estimated to be a 10% and a 15% 
respectively. 

On the contrary to the £ flavour, the 6 tagging 
efficiency depends on the number of jets in the 
event. The mean energy of the B hadrons in a 
three-jet event is smaller than that in a two-jet 
event. The same applies to their products and 
due to the multiple scattering, the error of their 
impact parameter is larger. Consequently, the 
probability of having all particles compatible with 
being produced in the IP is also larger. 

The background of the 6 sample comes mainly 
from those £ events having either not identified 
decays or 7 conversions in the innermost layers of 
the detector. The large impact parameter of their 
products pulls V towards low values. 

2.4. Jet rates 

The observed three-jet cross section of each 
tagged sample, q = b and £, was computed as: 

^3q {yc) — 3j-obslycj/-L tot-obs 

The measured value of our observable can be 
computed simply with the quotient of the recon- 
structed three-jet rates: 

This raw value of the observable must be con- 
verted into a parton level one which may be com- 
pared with the NLO calculations [0, || . The cor- 
rection method accounts for the detector effects, 
biases introduced in the flavour tagging plus the 
hadronization effects. 

The contribution of each flavour, R\^, to the 
observed three-jet cross section is given by: 

pobs b rjb I c pc I I r)l 

— • ^3q + ' -^3^ -|- • It^g 

where corresponded with the flavour contents 
for i = b, c, ^ of both tagged samples. These 
factors were extracted from the simulation. 

The reconstruction and parton level three-jet 
rates of each flavour and sample {R\„ and R^^'^ 



respectively) are related by: 

RUvc) = n^ivc) ■ ■ RTivc) 

where /g^ corresponds with the factor correct- 
ing both: detector acceptance and tagging effects. 
These factors can be deduced from the modelling 
of the DELPHI response to the hadronic events. 
Of course, data sets corresponding to different 
years of the detector's operation must be cor- 
rected independently. The hadronization factors, 
gsi, can be estimated by comparing the hadron 
and parton level distribution of the three-jet rate. 

The impact of the c quark mass in our observ- 
able can be safely neglected. As mass effects are 
proportional to m^, the effect of c quarks com- 
pared to &'s is {rric/mb)^ ^ 1/10 times a 10% 
factor, due to the contents on c flavour of the 
tagged samples. Hence, the influence of the c 
mass would be roughly an 1% that of the 6, i.e. a 
repercussion below the 0.5°/oo in both, numerator 
and denominator of the observable. Furthermore, 
the ratio reduces the net effect to negligible lev- 
els. According to this, the parton level jet rates 
of the c quarks are taken equal to those of the 
light quarks: R^^'^ = R^g" ■ Now, the measured 
jet rates can be expressed as: 

R°t\y,) = A,{y,) • R^^^y,) + B,{y,) ■ Rlf (y,) 
RfriVc) = A,{y,) ■ i?P-(ye) + B,{y,) • RlY{y,) 

where Aq and Bq are a redefinition of the orig- 
inal set of parameters: c*, /J^ and g^i- This 
parametrization allows expressing the corrected 
observable as: 

Tjpar 7, rj 7j6£-obs 

T^M, N _ -^3^ ^ Bb- Bi- R^ 

^3 yyc) ppar . nfc£-obs a 

The total correction applied to the raw R^ was 
about the 10%, the bulk of which corresponded to 
the detector plus tagging effects and an 1% to 
the hadronization. The corrected R^ distribution 
is shown in figure ^. The experimental points are 
seen to lay well below 1. So, effectively, h quarks 
radiate less gluons than light quarks do. In or- 
der to quantify this statement, the measurement 
must be based on a single point, since all data 
points are highly correlated. Large values of yc 
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Figure 1. Experimental values of i?3^(yc) compared 
with different QCD predictions at LO and NLO for 
two different mass hypothesis m;, = 3, 5 GeV/c'^. The 
experimental points are displayed only with their sta- 
tistical error. 

must be avoided in order to keep low the statisti- 
cal error. Also low values of yc should be eluded, 
thus limiting the effects of the higher order terms. 
Then, y^ = 0.02 is chosen and the result is: 

i?^^(0.02) = 0.971 ± 0.005(stat.) ± 0.007(frag.) 

2.5. Error estimation 

The correction of the measured RS^ induced 
some systematic uncertainties. The source of sys- 
tematics considered are the following: fragmenta- 
tion and simulation. 

The fragmentation model uncertainty has been 
evaluated by the comparison of two models: 
string fragmentation (JETSET [ |l||) and clus- 
ter decay (HERWIG [0). Huge samples were 
used (more than 10^ simulated events per gen- 
erator and flavour). The adopted value of the 
i?3^ was the average of those extracted correcting 
with both models, which differed by about 1%. 
The fragmentation model uncertainty was taken 
to be half of difference between them. 

The effect of the main fragmentation parame- 
ters (Qo, CTg, eb, a and h in JETSET) was stud- 
ied also with massive statistics (10^ events per 
flavour and parameter). Each of the above pa- 
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Figure 2. Experimental result of mb{Mz). The DEL- 
PHI point at Mz is displayed only with its statistical 
error. It is seen to be compatible with the running of 
nih as predicted by QCD. 



rameters took the optimum value found in the 
DELPHI tuning [ |l^ and varied by ±2a. For 
He > 0.01, none of the uncertainties due to the 
individual parameters exceeded the 2°/ao . The 
total error was obtained by adding quadratically 
the individual ones, and neglecting the correla- 
tion. It represented roughly a 3°/oo ■ 

The defects of the DELPHI detector modelling 
when trying to reproduce the real data were also 
considered in the analysis and regarded as an ad- 
ditional source of error. The main contribution to 
this uncertainty came from the limited statistics 
of the simulated events (4 • 10^) with fuU DEL- 
PHI simulation. Thus, the limited knowledge on 
the purity /background factors, , of the tagged 
samples entered as the simulation error. This er- 
ror has a strong dependence on yc as it is linked 
to the statistics of the three-jet simulated sam- 
ple. But, for yc = 0.02, it represents a mere 3°/oo 
(already included in the statistical error). 



3. The running mass of the b quark 

Using the DELPHI result of i?|^(yc) and the 
recent NLO calculations for exactly the same ob- 
servable [ , the value of the running mass of the 



b quark at the AIz scale is found to be: 

mb{Mz) = 2.67 ± 0.50 GeV/c^ 

where the breakdown of the error is listed in figure 
||. The theoretical uncertainty has two sources. 
The first one is the error due to the QCD scale. 
It has been estimated by varying the scale, /it, in 
the range: 0.5 < fi/Mz < 2, and represents 0.10 
GeV/c^ The second one is due to the mass ambi- 
guity. It accounts for the uncertainty of express- 
ing the mass dependence of R^^ by either using 
directly the running mass, rrib, or using the pole 
mass, Mb, and then run up to Mz using the RGE. 
Both practices are licit, but differ on the treat- 
ment of the higher order terms (with truncated 
or resummed expressions). The contribution to 
the error is 0.25 GeV/c^. 

The running of rub is verified when comparing 
its value at the Mz with that at the T scale. The 
net change in the b running mass between both 
scales is (errors added up quadratically): 

mbiMr/2) - mb{Mz) = 1.49 ± 0.52 GeV/c^ 

This result represents the first experimental ev- 
idence of the running property of the b quark 
mass, in particular, and of the any fcrmion mass, 
in general. 
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